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PCA+CHMM 在 设备 性 能 退化 状态 识别 中 的 应 用 研究 
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摘 要 : 为 了 准确 识别 机 械 设备 当前 所 处 的 退化 状态 ,研究 了 一 种 基于 PCA( 主 成 分 分 析 ) 和 CHMM (连续 型 隐 马 尔 
可 夫 模 型 ) 结 合 的 性 能 退化 状态 识别 方法 。 首 先 提取 设备 振动 信号 全 寿命 周期 的 时 域 、 频 域 、 时 频 域 的 特征 ， 经 过 初 
步 得 选 后 组 成 新 的 特征 集 ， 使 用 PCA 方法 对 其 进行 降 维 处 理 ; 然后 利用 降 维 后 的 数据 ， 训 练 一 个 全 寿命 周期 CHMM 
用 来 确定 退化 状态 数目 ， 再 针对 每 个 退化 状态 训练 一 个 CHMM,， 通过 比较 观测 序列 处 于 各 个 模型 下 的 似 然 概 率 值 判断 
设备 当前 所 处 的 退化 状态 ; 最 后 通过 实验 对 比 了 PCA+CHMM 和 PCA+SVM、PCA+KNN、PCA+CART 方法 的 各 退化 
状态 识别 准确 率 ， 结 果 表 明 PCA+CHMM 的 平均 识别 准确 率 最 高 、 识 别 效果 较 好 ， 适 用 于 设备 退化 状态 的 识别 。 
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Application of PCA + CHMM in equipment performance degradation state recognition 


Zhong Fei Ning Qian'?, Zhou Xinzhi!?*, Zhao Chengping! 
(1. College of Electronics & Information Engineering Sichuan University, Chengdu 610065, China; 2. Science & Technology 
on Electronic Information Control Laboratory, Chengdu 610036) 


Abstract: In order to accurately identify the degradation state of mechanical equipment, this paper researched a recognition 
method of performance degradation state based on PCA (principal component analysis) and CHMM (continuous hidden markov 
model) . Firstly, extracted the vibration signal’s time domain, frequency domain and time-frequency domain features in full life 
cycle, then constructed a new feature set by screening the features, then performed PCA dimensionality reduction for this set; 
Secondly, trained a full life cycle CHMM to determine the number of degraded states by using of reduced dimension feature 
data, And then training a CHMM for each degraded state, judging the degradation state of the device by comparing the likelihood 
probability of the observation sequence under each model; Finally, the accuracy of PCA + CHMM and PCA + CHM, PCA+ 
KNN and PCA+CART methods that to identify each degraded state are compared. The results show that the average recognition 
accuracy of PCA + CHMM is the highest, the recognition effect is good, and it is suitable for the identification of device degraded 
state. 

Key Words: performance degradation; principal component analysis; continuous hidden markov mode!l; feature dimensionality 
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设备 性 能 退化 状态 识别 的 主要 流程 包括 原始 数据 的 特征 提 
取 和 识别 模型 的 选取 和 建立 。 在 特征 提取 阶段 ， 由 于 单一 性 能 

IMS(intelligent maintenance systems) 概 念 由 美国 威斯康星 。 特征 参数 包含 信息 量 少 、 抗 干扰 性 差 、 故 障 敏感 度 低 、 信 息 宛 
大 学 Leel! 教 授 率先 提出 , 其 中 , 设备 性 能 退化 过 程 的 评估 与 预 。 余 等 外 缺点 ， 一 般 会 同时 提取 多 个 特征 参数 进行 设备 的 性 能 衰 
测 是 IMS 的 核心 技术 。 正确 识别 设备 当前 所 处 的 性 能 衰退 状态 。 退 评 佑 。 特 征 参 数 维 度 很 高 时 不 可 避免 地 会 包含 一 些 见 余 信 息 
是 进行 预测 性 维护 的 前 提 条 件 ， 将 性 能 退化 识别 技术 应 用 到 实 。 和 噪声 ， 使 得 模型 训练 容易 出 现 过 拟 合 的 现象 ， 所 以 必须 采 上 
际 的 工业 生产 中 去 ， 可 以 极 大 程度 地 减少 设备 的 维修 时 间 和 避  ” 特征 降 维 方法 对 特征 数据 进行 降 维 。 特 征 降 维 可 通过 特征 选择 
免 设 备 故障 引起 的 灾难 性 事故 的 发 生 ， 对 降低 维护 成 本 和 提高 (feature selection) 或 者 特征 抽取 (feature extraction) 实 现 申 ， 特 征 
生产 效率 具有 十 分 重要 的 意义 。 选择 只 是 简单 地 从 原 特 征集 里 面 选取 子 特 征集 ， 而 特征 抽取 则 
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是 对 原 特征 集 的 某 种 映射 变换 。PCA 就 是 一 种 常用 的 特征 


算法 ， 


已 能 在 保证 数据 信 


取 
息 丢 失 最 少 的 原则 下 用 更 少 的 不 相关 
关 


、 
各 


综合 变量 蔡 代 原本 较 多 的 变量 ; HMM ( 隐 马 尔 可 夫 模 型 ) 是 


于 时 序 的 概率 模型 ， 在 语音 识别 、 自 然 语 言 处 理 、 模 式 识别 等 


领域 有 着 广泛 的 应 | 
式 分 类 与 识别 中 存在 


]， 考 虑 到 设备 的 故障 诊断 和 语音 识别 在 模 


定 


的 相通 性 ， 因 此 HMM 可 以 用 于 设备 


的 性 能 退化 过 程 建 模 外 ,HMM 定义 的 观测 数据 类 型 为 离散 值 ， 


当 观 测 值 为 连续 型 数 扫 


届时 ， 先 要 对 特征 集 进 行 矢量 量化 才能 进 


型 的 识别 精度 。 CHMM 使 用 混合 高 共 
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的 和 矩阵， 计算 各 主 分 量 的 贡献 率 (贡献 率 = 单个 特征 值 /特征 值 
总 和 ) ， 按 需 选 取 其 中 部 分 主 分 量 。 
1.2 CHMM 定义 

HMM ( 隐 马 尔 可 夫 模 型 ) 最 初 由 Baum 等 人 提出 ， 它 是 一 
个 关于 时 序 的 概率 模型 ， 具 有 很 强 的 时 序 模式 分 类 能 力 。HMM 
是 双重 随机 过 程 : 一 个 是 隐 含 马尔 可 夫 链 的 状态 转移 过 程 ， 一 
个 是 描述 隐 含 状态 与 观测 序列 之 间 关 系 的 观测 过 程 。 不 同 于 
HMM，CHMM 的 观测 序列 为 连续 值 ， 观 测 概率 矩阵 由 混合 高 


密度 函数 来 描述 连续 型 观 


载 设备 上 
换 和 隐 马 尔 可 夫 模 型 
取 设 备 振 
对 轴承 怕 


的 分 析 实 验 订 


1 基本 理论 


1.1 PCA 介绍 


PCA (principal component analysis) 是 一 种 较为 常用 的 特征 


寺 

行 HMM 的 建 模 ， 这 势必 导致 特征 信息 的 丢失 ， 从 而 影响 到 模 
折 4 
由 


测 值 的 分 布 ， 可 以 直接 利用 连续 观 涡 
矢量 量化 带 来 的 信息 丢失 。 文 献 [5] 采 用 PCA 和 DHMM (离散 
型 隐 马 尔 可 夫 模 型 ) 结合 的 方式 应 用 于 发 动机 的 故障 诊断 ; 文 
献 四 引入 状态 条 件 概 率 矢 量 对 HMM 进行 不 确定 性 改进 后 对 机 
的 健康 状态 进行 评 
4 结合 
动 信号 的 小 波 包 粹 作为 特征 向 量 ， 结 合 高 斯 混合 模型 
E 能 退化 过 程 进行 评估 ;文献 [9] 成 功 将 EEMD 和 SVM 
结合 起 来 应 用 到 滚动 轴承 的 退化 状态 识别 ;文献 [10] 则 采用 隶 
属 度 函数 和 CHMM 对 轴承 的 性 能 退化 过 程 建 模 ， 取 得 了 较 好 
的 评估 效果 。 本 文 研究 了 一 种 PCA 特征 降 维和 CHMM 结合 的 
设备 性 能 退化 状态 识别 方法 ， 并 通过 一 个 全 寿命 周期 轴承 数据 
E 明 了 此 方法 拥有 较 好 的 识别 效果 。 


降 维 方法 ， 由 Pearson[ 于 1901 年 提出 。 主 成 分 分 析 的 基本 原 
里 是 通过 一 个 正 交 化 线性 


分 量 不 相关 的 新 随机 向 


蔡 原 随机 变量 ， 


三 


量 


值 进行 模 型 训练 ， 避 免 J 


十 ; 文献 [7] 提 出 了 一 种 基于 小 波 包 变 
的 轴承 性 能 退化 评估 方法 ; 文献 [8] 提 


斯 密度 函数 组 成 。 
CHMM 可 以 由 下 面 四 个 参数 描述 [1: 
引 隐 含 的 状态 个 数 N， 可 将 隐 含 状态 分 别 记 为 3 、S> 、 
S3 ee SN， 时 刻 t 的 状态 记 为 9 ,那么 有 4,€(51,52,53.…… SN) ; 
b) 初 始 时 刻 1=0 时 的 各 状态 分 布 概率 
X= (N,N A3...... ZN) ， 简 记 为 T=P(qgg=5;),，1<i<N; 
9] 隐 含 状态 之 间 互 相 转移 的 概率 矩阵 A={aijNxN ，4 就 是 
状态 5; 转移 到 8 的 条 件 概率 值 ， 记 为 中 =PP(4qiw =Sj|4=5)， 
1<i,j<N; 


变换 将 分 量 相 关 的 原 随机 向 量 转 换 为 
， 新 随机 向 量 的 分 量 按 方差 贡献 率 的 


主 分 量 


区 其 中 综合 贡献 率 超过 85% 的 部 分 主 分 量 就 可 以 代 
达到 数据 降 维 的 目的 。PCA 在 最 大 程度 地 保留 


、 第 二 主 分 量 、 第 三 主 分 量 ， 依 次 类 


原 随机 向 量 信息 的 前 提 下 , 消除 了 特征 参数 之 间 的 相关 性 .PCA 


ES 


的 计算 步骤 描述 如 下 : 
a) 将 原始 特征 集 的 每 个 特征 作为 列 向 量 ， 构 造 N*M 的 样 


d 观 测 值 概率 矩阵 B=fbi(Oojx ， 其 中 ， 


M 
bj(0D) = cmN(orspimU im) ，1< j<N，M 代表 混合 高 斯 元 


m=1 


的 个 数 ，Cjn 是 混合 高 斯 元 的 权重 ，Kim 是 均值 向 量 ,Uj 是 协 


M 
方差 矩阵 ， 同时 满足 : > cm>0 和 | b;(Wadx=1; 


m=1 
综 上 所 述 ，CHMM 可 以 表示 成 4=(z,N,A,cjns HjmsU jm) 。 
CHMM 包括 三 个 基本 问题 : 评估 问题 、 解 码 问题 和 学 习 问 
题 ， 相 对 应 的 算法 分 别称 为 前 向 一 一 后 向 算法 、Viterbi 算法 和 
Baum-Welch 算法 。 算 法 相关 的 详细 推理 过 程 同样 可 参考 文献 
[9]。 


2 PCA+CHMM 退化 状态 识别 方法 


设备 健康 状态 从 正常 到 失效 会 经 历 一 系列 的 退化 状态 ， 对 
退化 状态 的 识别 实质 上 是 模式 识别 问题 ， 常 见 的 模式 识别 算法 
有 KNN、 朴 素 贝 叶 斯 、 决 策 树 和 SVM 等 。 如 果 能 正确 识别 出 
设备 当前 时 刻 的 运行 状态 ， 就 可 以 在 故障 发 生 之 前 采取 相应 的 
预防 或 者 修复 措施 , 从 而 避免 重大 故障 的 发 生 , 减少 维护 成 本 ， 


本 矩阵 X ，NN 为 特征 长 度 ，M 为 特征 集 的 特征 个 数 ; 


减 去 对 应 的 均值 得 到 新 


各 分 量 再 除 以 N-1); 
9 求解 X 的 M*M 协 方差 矩阵 C ; 


d) 求 解 协 方差 外 


入 


b) 对 X 进行 中 心 化 处 理 ( 和 矩阵 中 心 化 就 是 先 将 每 个 列 向 量 
E 阵 ,然后 新 矩阵 和 新 矩阵 的 转 置 相 乘 ， 


E 阵 C 的 特征 值 和 特征 向 量 ; 


将 第 中 步 得 到 的 特征 值 按 从 大 到 小 的 顺序 排列 并 构造 成 一 


个 列 向 量 ， 依 此 顺序 将 对 应 的 特征 向 量 作为 列 向 量 组 合成 特征 
向 量 矩阵 ; 


e) 匀 


E 阵 X 和 特征 向 


地 


利用 PCA 和 CHMM 结合 的 方法 就 可 以 实现 对 设备 性 能 状态 的 
准确 识别 。 进 行 设备 退化 状态 识别 之 前 需要 采集 设备 的 运行 数 
据 ， 传 感 器 获取 的 原始 信号 数据 量 大 且 包 含 噪声 ， 所 以 有 必要 
对 原始 数据 进行 预 处 理 。 

本 文 提出 的 退化 状态 识别 方法 流程 描述 如 下 : 

a) 对 设备 原始 数据 进行 预 处 理 ， 为 保证 提取 特征 的 全 面 性 
和 有 效 性 ， 分 别 从 时 域 、 频 域 和 时 频 域 三 个 方面 进行 提取 。 时 
域 特征 分 为 有 量 纲 和 无 量 纲 两 种 ， 有 量 纲 特征 包括 均值 、 均 方 


和 矩阵 相 乘 就 得 到 PCA 的 主 分 量 组 成 


根 值 、 方 根 幅 值 、 绝 对 平均 值 、 牌 度 、 峭 度 、 方 差 、 最 大 值 、 
最 小 值 、 峰 峰值 10 种 ， 无 量 纲 特征 包括 波形 指标 、 峰 值 指标 、 
脉冲 指标 、 裕 度 指 标 ， 峭 度 指 标 和 偏 斜 度 指标 6 种 ， 提 取 频 域 
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特征 时 先 对 原始 信号 做 快速 傅 里 叶 变换 ， 再 / 


频谱 和 各 频率 分 


量 值 构造 共 13 种 特征 ， 顺 序 编号 为 频率 特征 1-13; 时 频 域 特 
征 是 原始 数据 的 非 线 性 特征 ，EMD (经 验 模 态 分 解法 ) 能 够 提 
取出 原 信 号 中 的 微弱 特征 ， 利 用 EMD 将 原始 信号 分 解 为 6 个 
IMF 〈 本 征 模 态 函数 ) 分 量 ， 各 IMEF 分 量 的 能 量 值 组 成 6 种 时 
频 域 特征 。 以 上 特征 共计 35 种 组 成 特征 集 , 各 特征 参数 的 定义 
式 可 参考 文献 3]; 

b)PCA 特征 降 维 ， 从 上 一 步 获取 的 特征 集中 筛选 出 能 够 明 
显 体现 设备 退化 趋势 的 特征 组 成 新 的 特征 集 ， 按 照 1.1 小 节 描 


述 的 步骤 对 新 的 特征 集 进行 PCA 降 维 ， 
贡献 率 ， 选 取 综 合 贡 献 率 超过 一 定 阔 


样本 数据 ; 
oO) 确定 设备 退化 状态 数 

ee spr 

入 置 为 [4,10]， 


再 利用 


Viterbi 


为 最 终 的 退化 状态 数目 ; 


d) 针 对 不 同 的 退化 状态 训练 不 同 的 CHMM， 共 训 
计算 当前 观测 序列 在 这 N 个 模型 下 的 似 然 概率 值 ， 


通过 


， 有 具体 做 法 是 用 
CHMM， 全 周 


并 分 别 计算 各 主 分 量 的 
值 的 主 分 量 组 合成 最 后 的 


上 一 步 获 取 的 所 


解码 ， 选 取 退 化 


期 CHMM 的 隐 
算法 进行 


状态 划分 误差 小 且 符合 设备 退化 状态 时 序 转移 规律 的 状态 数 作 


练 N 个 ， 
似 然 概 


率 值 大 的 模型 就 是 当前 设备 的 退化 ; 


提取 时 域 、 频 域 、 时 频 域 的 35 种 特征 


筛选 有 效 特 征 ， 进 行 PCA 特 征 降 维 


训练 全 寿命 周期 CHMM 


池 注 讲演 民 遍 字 联 了 流 毛 T41179 人 


当前 退化 状态 


图 1 


PCA+CHMM 设备 性 能 


3 ”实验 分 析 


数据 说 明 
实验 采用 美国 辛辛那提 大 学 智 
Bearing 1 的 全 寿命 周 


3.1 


局 


采样 后 轴承 失效 。 如 图 2 所 示 。 


退化 状态 识 


期 加 速 疲劳 振 
2115， 实 验 时 转速 2000 r/min ， 径 向 载荷 26.6 kN ， 
20480 Hz， 每 10 分 钟 采 样 一 次 , 每 次 采样 持续 1 秒 ， 
原始 数据 


大 态 。 


运行 监测 数据 


训练 退化 状态 1 的 CHMM 


训练 退化 状态 2 的 CHMM 


训练 退化 状态 N 的 CHMM 


比较 似 然 概 率 值 Jog(P(O|4)) 


别 方法 流程 图 
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第 984 次 
不 能 清晰 地 反映 轴 


承 的 性 能 退化 趋势 。 
3.2 ”特征 提取 和 PCA 降 维 


从 时 域 、 频 域 、 时 频 域 提取 原始 信号 


的 35 种 特征 ,筛选 出 


对 退化 状态 变化 敏感 的 特 4 


正 组 成 特 和 
值 随 着 时 间 变 化 ， 具 有 一 定 的 趋势 规 


E 集 。 例 如 ， 时 域 的 均 方 根 
律 ， 而 至 度 特征 则 在 设备 


、 
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失效 前 几乎 看 不 出 有 什么 变化 ， 显 然 在 特征 
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图 3 时 域 特征 选取 对 比 图 
同样 的 方法 ， 对 频 域 和 时 频 域 的 特征 进行 筛选 。 最 后 选 
出 了 均值 、 均 方 根 值 、 方 根 幅 值 、 绝 对 平均 值 、 最 大 值 、 最 小 
值 、 峰 峰值 、 波形 指 标 、 偏 斜 度 指 标 、 峭 度 指 标 、 频 域 特征 1、 
频 域 特征 2、 频 域 特征 6、 频 域 特征 10、 频 域 特征 13、IMF1 能 
量 值 、IMF3 能 量 值 、IMF4 能 量 值 、IMF5 能 量 值 、IMF6 能 量 
值 共 20 种 特征 , 此 时 特征 维 数 仍然 较 高 ， 如果 将 其 直接 作为 样 
本 数据 ， 一 方面 数据 量 大 ， 处 理 起 来 不 方便 ， 男 一 方面 很 多 特 
征 的 趋势 非常 相似 ， 存 在 信息 元 余 
0.05 0.9 
IMF4 能 量 值 IMF6 能 量 值 
0.045 | 0.8|- 
0.04|- | 
0.035 
0.6|- 
0.03|- 
0.5|- 
0.025 
0.4 上 - 
0.02|- 
0.3|- 
0.015 | 
0.01 上 | 0.2 | 
0.005 0.1 机 
上 | "| 
0 200 400 600 800 1000 0 200 400 600 800 1000 
特征 时 序 点 特征 时 序 点 
到 4 ”时 频 域 特征 存在 相似 性 (元 余 ) 


对 这 20 种 特 生 
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正 进 行 PCA 降 维 ， 


得 到 第 
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0.02%， 前 两 个 主 分 量 综合 贡献 率 达 99.96%， 


所 以 选取 第 一 和 


第 二 主 分 量 作 为 最 终 的 样本 数据 特征 
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图 5 PCA 降 维 后 的 特征 


3.3 CHMM 建 模 与 测试 


上 一 小 节 PCA 降 维 得 到 的 全 周期 样本 数据 作为 训练 序列 ， 
分 别 设置 隐 舍 状态 数 为 4 至 10， 利 用 Baum-Welch 算法 训练 得 
到 不 同 的 全 周期 CHMM， 再 分 别 采 用 Viterbi 算法 对 样本 序列 
进行 解码 ， 通 过 对 比 退 化 状态 划分 效果 来 确定 退化 状态 数目 。 
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图 7 状态 数 为 5-10 时 的 划分 效果 


如 图 ， 


状态 数 设置 为 5、6、 


所 以 将 设备 退化 状态 数 设 置 为 4， 划 分 
将 全 周期 样 


7 时 ， 最 后 一 个 状态 的 数据 点 
只 有 1 个 ， 显 然 不 符合 实际 ;状态 数 设置 为 8,、9、10 时 
分 结果 出 现 交 错 现象 ， 也 不 符合 设备 性 能 的 退化 规 得 
态 数 设置 为 4 时 的 划分 结果 较为 符合 性 能 退 
区 间 见 表 1。 
本 数据 分 成 4 个 状态 区 间 ， 每 个 状态 区 间 的 样 
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本 数据 均 按 7:3 的 比例 划分 训练 序列 和 测试 序列 ，4 个 训练 序 
列 训练 得 到 4 个 不 同 退化 状态 下 的 CHMM ， 


前 观测 


序列 在 每 个 退化 状态 ae 民 化 状 
态 的 识别 。 
表 1 设备 性 能 退化 状态 区 间 划 分 及 CHMM 数据 划分 
状态 正常 状态 ” 轻 度 退化 ”中 度 退 化 ”重度 退化 
区 间 [1,700] [701,879] [880,957] [958,984] 
训练 序列 [81,570] [731,855] [891,944] [962,980] 
测试 序列 210 组 54 组 23 组 8 组 


如 表 1 所 示 ， 分 别 将 各 个 状态 的 测试 序列 输入 到 4 个 模型 


中 ， 


得 到 4 种 似 然 概 率 输 


出 值 的 对 比如 图 
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8~11 所 示 。 
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图 11 


E 度 退化 测试 集 的 似 然 概率 值 比较 


Lm 


同样 地 ， 将 全 寿命 周期 样本 数据 《〈 特 征集) 按 7:3 划分 为 
训练 集 和 测试 集 ， 分 别 训练 KNN、CART、SVM (多 分 类 ) 模 
型 ， 对 轴承 的 性 能 退化 状态 进行 识别 , CHMM、KNN、CART 
和 SVM 识别 准确 率 对 比如 表 2 所 示 。 


表 2 轴承 退化 状态 识别 效果 对 比 
- 正常 轻 度 中 度 重度 平均 
状态 、 、 网 

态 退化 退化 退化 准确 率 

PCA+CHMM 99.52% 96.30% 100% 100% 98.21% 
PCA+SVM 99.52% 100% 86.96% 75% 90.37% 
PCA+CART 97.62% 92.59% 95.65% $50% 83.97% 
PCA+KNN 98.10% 96.30% 100% 75% 92.35% 


表 2 中 ，PCA+CHMM 在 四 种 方法 里 退化 状态 识别 的 平均 
准确 率 最 高 , 且 每 种 状态 的 识别 准确 率 均 超过 96%, 比较 均衡 ; 
PCA+CART 的 平均 识别 准确 率 最 低 , 它 对 重度 退化 状态 的 识别 
效果 非常 差 .， PCA+SVM 和 PCA+KNN 的 平均 识别 准确 率 均 在 
90% 以 上 ， 但 是 对 各 退化 状态 的 识别 准确 率 却 不 是 很 均衡 。 不 
同 于 KNN、SVM 和 CART 等 模式 识别 算法 , CHMM 可 以 为 每 
个 退化 状态 训练 一 个 模型 ， 具 有 较 强 的 时 序 模式 识别 能 力 ， 恰 
好 符合 设备 退化 状态 在 时 间 上 分 布 不 均衡 的 情形 。PCA 和 
CHMM 结合 的 方法 在 设备 性 能 退化 状态 识别 上 能 够 取得 了 较 
高 的 识别 准确 率 ， 实 验证 明 此 方法 可 行 且 有 效 。 


4 ”结束 语 


本 文 提出 PCA 和 CHMM 结合 的 设备 性 能 退化 状态 识别 方 
法 ，PCA 起 到 了 去 除 信息 见 余 、 降 低 数据 维度 的 作用 , CHMM 
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不 仅 具 有 很 强 的 时 序 模式 识别 能 力 ， 而 且 避 免 了 矢量 量化 带 来 
的 信息 损失 。 利 用 美国 辛辛那提 大 学 提供 的 轴承 数据 ， 通 过 实 
验 对 比 知道 ， 此 方法 适用 于 设备 退化 性 能 状态 的 识别 ， 具 有 较 
高 的 识别 准确 率 。 因 为 此 方法 只 能 识别 当前 设备 所 处 的 退化 状 
态 ， 未 来 的 工作 可 以 从 回归 算法 和 CHMM 结合 的 方向 着 手 ， 
实现 对 设备 性 能 退化 状态 的 预测 。 
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